Hawking radiation in sonic black holes 
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I present a microscopic description of Hawking radiation in sonic black holes. A one-dimensional 
Fermi-degenerate liquid squeezed by a smooth barrier forms a transonic flow, a sonic analogue of 
a black hole. The quantum treatment of the non-interacting case establishes a close relationship 
between the Hawking radiation and quantum tunnelling through the barrier. Quasi-particle exci- 
tations appear at the barrier and are then radiated with a thermal distribution in exact agreement 
with Hawking's formula. The signature of the radiation can be found in the dynamic structure fac- 
tor, which can be measured in a scattering experiment. The possibility for experimental verification 
of this new transport phenomenon for ultra-cold atoms is discussed. 
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Introduction. — One of the most fascinating effects in 
astrophysics is Hawking's prediction of black hole radia- 
tion Hawking realized that particles such as neutri- 
nos and photons could be created spontaneously from any 
black hole and emitted at just the rate that one would 
expect if the black hole is a body with a temperature 
hK,/2-KkB proportional to its surface gravity k. Although 
Hawking's prediction combines gravitation and quantum 
mechanics to produce thermodynamics, it was rapidly 
appreciated that the radiation predicted by Hawking is 
more primitive and fundamental. A beautiful example 
of this was pointed out by Unruh 0|: The propagation 
of sound in a fluid/gas turning supersonic, such as in 
the nozzle of a rocket engine, is similar to the propaga- 
tion of a scalar field close to a black hole. By adapting 
Hawking's derivation, Unruh predicted the analogue of 
Hawking radiation for sound waves in a moving fluid. 
This sonic analogue could provide a way to experimen- 
tally test Hawking's effect. 

Although Unruh's prediction is not restricted a pri- 
ori to a specific fluid, superfluids immediately were con- 
sidered as the natural candidates. Superfluid helium II 
has been the first system to be proposed for the inves- 
tigation of the Hawking effect 0. More recently atomic 
Bose-Einstein condensates have been considered, as they 
may offer a better chance of testing Hawking radiation 
El El Transonic flows could be also realized in one- 
dimensional (ID) systems 0, 0] and available ID quan- 
tum systems, including ultra-cold atoms in optical po- 
tentials and wave-guides |nj > could be adapted to model 
a sonic black hole analogue and to investigate the elusive 
Hawking effect. 

The Hawking effect has not been considered so far in 
systems of non- interacting fermions. It is the aim of 
this Letter to show that indeed a ID Fermi-degenerate 
non-interacting gas that scatters against a very smooth 



potential barrier provides a clear and straightforward 
quantum mechanical microscopic description of Hawk- 
ing radiation. Sonic Hawking radiation is formed due 
to quantum tunnelling through the top of the potential 
barrier. Particle-hole excitations are created at the bar- 
rier and are then radiated with a momentum distribution 
characterized by a temperature in exact agreement with 
Hawking-Unruh's formula 0, 0, 0] ■ Although the evolu- 
tion of the fluid's flow is unitary, when the momentum 
distribution is measured locally, for instance through the 
measurement of the dynamic structure factor, it is not 
distinguishable from a thermal (and incoherent) distri- 
bution. The general case of ID quantum fluids as well as 
the possible relevance for experiments is also discussed. 

Part I. Microscopic description. — Consider a Fermi- 
degenerate gas of particles moving in a sufficiently nar- 
row and smooth constriction such that their motion can 
be considered ID. Their stationary single-particle wave- 
functions can be written in the adiabatic approxima- 
tion [ill as W(x,y,z) = ip(x)£ . x (y, z) where £o- tX (y,z) is 
the transverse ground-state wave-function with eigenen- 
ergy eo(x). ip(x) satisfies a ID Schrodinger equation 
[— (H 2 /2m)d 2 . + V ex t(x)]tp(x) = eijj(x) where V cxt is an ef- 
fective potential that includes also the zero point energy 
€q(x) [lJI- I assume V^xt very smooth and essentially 
non-zero only near its maximum V max located at x — 0. 

Consider a flow of particles coming from the left and 
colliding against the potential V ox t. Their particle's wave- 
functions are given by scattering wave functions asymp- 
totically defined by ip e oc [e zkx + r(e)e~ lkx ]Q(-x) + 
t(e)e tkx Q(x), where e = h 2 k 2 /2m is the energy eigenvalue 
and r(e) and t(e) are the reflection and transmission am- 
plitudes, respectively. The evolution of the particle flow 
is unitary and the many-body wave function of the flow 
of particles is represented in a second quantization by 
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where a| is the creation operator of V'e's and |— ) is the 
vaccum of particles. In this state all ^ e 's are occupied 
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up to an energy e max , which is assumed higher than 
Knax' The non- interacting case is an ideal example of 
a fluid turning supersonic (see below) where the micro- 
scopic wave function of the flow JIJ is given from the 
beginning. 

Since the barrier potential is smooth the particle 
propagation is semi classical and in first approximation 
|t(e)| 2 w 0(e — V max ). Therefore the fluid maintains in 
a first approximation the sharp Fermi-degenerate char- 
acter of the local momentum (or velocity) distribution, 
which represents the "Vacuum" of the black hole sonic 
analogue. This allows us to describe the fluid dynamic 
with the standard hydrodynamic equations. The veloci- 
ties of the particles in the fluid are uniformly distributed 
within an interval (vl,vr). Figure f shows an example 
of velocity distribution. The right Fermi velocity vr = 
V "max — 214 x t (x) / m is the velocity of a particle mov- 
ing from the left to the right with initial velocity higher 
than the cla ssical thre shold velocity to go over the bar- 
rier w osc = y t /2V max /m. The left Fermi velocity vl is the 
velocity of a particle starting from the top of the barrier 
(located at x = 0), i.e. v L — sgn(ir) \/v% sc — 2V cxt {x)/m. 
The sgn(x) function accounts for the two possible differ- 
ent signs of the velocity that the particle can have falling 
either to the right or left of the barrier. 




FIG. 1: Semi classical velocity distribution as function of 
the coordinate x. Velocities are uniformly distributed in the 
area limited by Vr (top solid line) and vl (bottom solid 
line). Thinner solid lines correspond to different classical 
orbits of the particles of the fluid. Here e max = 2.5 and 
Voxt = exp(— x 2 ) in dimensionless units. The inset shows 
the Mach number v/c as function of position x. 



The flow velocity v = (vl + vr)/2 is the velocity of 
a local reference frame from where the fluid appears 
in equilibrium. Introducing the Fermi velocity vp = 
(— i>l + wr)/2 as the Fermi velocity in the co-moving 
frame, the density is given by n — mvp/izK. With these 
definitions the current — nv = mvvp/'K% is conserved 
[l3|. The Bernoulli equation v 2 /2 + h + V cx t/m = cost 
where h(n) = TT 2 h 2 n 2 /2m 2 is the enthalpy (the Fermi en- 
ergy divided by the particle mass) is satisfied with cost = 
v 2 sc /2 + z^ ax /2. The sound velocity (c 2 = ndh{n)/dn) 



is consistently given by the Fermi velocity c — vp. It is 
easy to verify that the above defined particles flow is a 
transonic flow, namely subsonic on the left of the barrier 
and supersonic on the right. The inset of Figure I shows 
the behaviour of the Mach number v/c in the ID channel. 

Even if the potential barrier is very smooth, its finite 
thickness near its top smears the semi classical Fermi 
distributions at vl(x). This is due to those particles col- 
liding against the barrier with energies comparable with 
Vmax- The reflection probability is given by (see also Q) 
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1 + cxp(27r(e - V max )/hw x ) 



(2) 



where oj x is the frequency of the inverted harmonic os- 
cillator obtained expanding the barrier potential V ex t = 
^muj 2 x 2 + o(x 2 ) around its maximum. Hawking's 



V max - 2 mUJ x- 

formula for a sonic black hole 



is given by 
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where the surface gravity k in Hawking temperature 
Kh/2nkB is replaced by d(v — c)/dx. For the transonic 
flow of fermions defined above d(v — c)/dx = d(vL)/dx = 
lj x and Hawking temperature 10 becomes 



huj x 

2vrfc F 



(4) 



Therefore the particles reflected from the barrier are dis- 
tributed according to the reflection probability J5J , which 
is identical to a thermal Fermi distribution 
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l + exp((e- M )/A:BT) 



(5) 



where the chemical potential /i = V max and the temper- 
ature is in exact agreement with Hawking temperature 

Consider a probe beam intersecting our ID flow in a 
region quite far from the barrier. The measurement of 
the differential scattering cross-sections with the probe 
beam allows obtaining the dynamic structure factor |15| , 
which can be used to investigate Hawking radiation. 
The dynamic structure factor is defined as S(q, to) — 
Zv f {f\p q \l){l\pl\f)5(u>-ui) where w, = and 
Ei are the excitation frequencies and the energies of the 
eigenstates \l), respectively. Positive (negative) frequen- 
cies in S(q,uj) are related to absorption (release) of en- 
ergy from the system. = £\ e }i(xi-Rp) e -(xi-Rp) /L p j g 
a density operator defined locally where L p ^> 1/q. Here 
L p and R p are the size and the position of the region 
where the probe beam and the system beam interact, 
respectively 

The dynamic structure factor of the flow is given in 
the subsonic region by 



S{q,u) = S^ R {q,u)0 [ - 
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and in the supersonic region by 



and /i/f, = e n 



e(-) (7) 



S^?,,(g, k>) is the 



where /x^ — v max an<a p.^ — c max . ^t.u 
dynamic structure factor of an equilibrium ID Fermi 
gas at temperature T and chemical potential /x given 
by SP )M {q,w) = ^ e ^/ fc B r /{Cosh(^/fc B T) + 



Cosh[(? 



^ + ^ - M)/^ B T]} where L cS - y n,^ p 
is an effective length of the probe-system interaction 
region [l6]. Its T = limit is given by S^^q^Lo) — 



e ( - u + ^1 + ) e 



2nK 2 q w I *~ 1 2m 

with k° F — ^/2mfi/h. 
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FIG. 2: Sonic Hawking radiation can be detected through the 
measurement of S(q,cj). The figure shows the w-dependence 
of S(q,u>) in the subsonic (left) and supersonic (right) re- 
gions and for a positive (top) and a negative (bottom) value 
of q. Here k B Tn = 0.15 V max , e max = 2.5 V max and \q\ = 
0.05 mv csc /h is chosen in the region of small q where thermal- 
like transitions (broadened peaks) dominate over the zero- 
temperature contributions (step-like peaks). S(q,u>) is in 
units of L e gm / 2ivh 3 q and u) is in units of Vm^/h. 



Due to the Pauli exclusion principle and to the one- 
dimensionality, which disconnects the Fermi- "surface" 
into two points (represented by ul and vr), S(q,u>) 
has a typical peaks structure, which contrasts the non- 
interacting two and three dimensional cases [13. The 
peaks are related to the dispersion relation u)(q) — 
—Vl \q\ 9(— q) + vr \q\ @{q) valid for small q. In the 
supersonic region Lo(q) is negative for q < meaning 
that small perturbations propagate always to the right. 
Figure 2 shows the typical behaviour of S(q,u>). The 
broad peaks in the bottom of Figure 2 (q < 0) corre- 
spond to the negative-g branch of the dispersion relation 
(oj(q) = —vl \q\). Also the broad peaks in the top of 
Figure 2 (q > 0) correspond to the same dispersion re- 
lation but with opposite sign of lo and q, but represent- 
ing these the inverse processes. In contrast the sharp 
peaks in the top of Figure 2 correspond to the dispersion 
relation Lo{q) — v^\q\, which is not broadened by the 



tunnelling through the barrier. The characteristic wave- 
vectors of Hawking radiation is q c ~ /cbTh/c ~ l/27r 2 ^n 
where l x — ^m/hu! x and n is the density. Indeed at 
very low frequencies fku <^ /sb?h and momentum S(q, lo) 
is dominated by Hawking radiation. The peak values 
at \u>\ = \v-^q\ tends to 1/4 of the corresponding zero- 
temperature value L e g m / litTc'q but with a half- width 
Auj = 1.76 qkBTn/mvL, which is much larger than the 
corresponding half- width Auj = hq 2 /2m of the zero- 
temperature contribution. 

The experimental observation of the widths of these 
peaks would provide a measure of Hawking tempera- 
ture. Degenerate Fermi atoms have been only recently 
obtained Lithium can have a degeneracy tempera- 
ture as high as 8 fiK, which fixes the scale of e max . A 
reservoir of degenerate Lithium could be connected to 
non trapped states with a ID channel realized by optical 
potentials or wave-guides 9] . The barrier potential could 
be created by a far-blue-detuned laser sheet propagating 
in a direction orthogonal to that of the flow propagation 
and tightly focused with a waist (1/e 2 ) of twice the laser 
wavelength (760 nm). Hawking temperature would be of 
order of Th ~ 200 nK and the characteristic wavelength 
of Hawking radiation A c of order of 40 fxxsx. 

The dynamic structure factor S(q, lo) in the subsonic 
region is equivalent to that of a ID channel connecting 
two reservoirs at temperatures T — on the left and 
T = Th on the right. This is the analogue of Hawking 
radiation from a black hole into the surrounding zero- 
temperature space. To distinguish between this pure 
quantum state and a mixed thermal state it would re- 
quire the measurement of correlations such as for instance 
those between opposite sides of the barrier. Correlated 
pairs of quasi-particles (particle- hole), the "real" parti- 
cles in an effective low-energy description, are created at 
the barrier and evolve in opposite directions. 

The interpretation of the Hawking radiation as related 
to a tunnel effect is in agreement with a recent derivation 
of the Hawking radiation as a tunnelling process of parti- 
cles in a dynamical geometry due to Parikh and Wilczek 
[l7| . A relationship with a generalized form of diffraction 
has been pointed out by Sanchez 0] . Moreover electric 
pair production in a strong static electric field is related 
to a similar tunnelling through an inverted harmonic os- 
cillator [Is[ . As quantum tunnelling is very sensitive to 
the details of the potential barrier, the distribution of ra- 
diated quasi-particles, related to the reflection coefficient 
of the barrier, can be different from J2J. In particular 
this suggest the possibility that some information on the 
supersonic part of the fluid (the analogue of the classi- 
cal forbidden region of a black hole) is carried by the 
radiation. 

Part II. ID quantum fluids. - Hawking's predic- 
tion is, in principle, valid also for strongly correlated 
ID systems, such as the Tonks and Lieb-Liniger gases 
[20| . which are currently subject of intense interest 01 
and electrons in clean short quantum wires PH2. us- 
ing standard hydrodynamics of ID flow, it is possible to 
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rewrite Hawking temperature (|3J in a more general form 
|2l| as T"h = r/ huj x /2TrkQ where n is defined by 




In general rj is only weakly dependent on the form of the 
function enthalpy h[n) and Hawking temperature is ba- 
sically proportional to a quantity that is related to the 
propagation of a single particle of the fluid in the vicinity 
of the horizon. In the case of a non-interacting degener- 
ate ID Fermi gas considered here, 77 is equal to 1 since the 
enthalpy is proportional to the square of the density. In 
the case of a ID Bose gas, 77 is equal to \/3/4 in the mean- 
field regime since the enthalpy is linear in the density 
and is equal to 1 in the Tonks-Girardeau limit [2(| since 
the enthalpy is identical to that of a Fermi-degenerate 
gas in that limit. Beside the importance of testing some 
aspects of black hole evaporation, sonic Hawking radia- 
tion would be a remarkable transport phenomena in these 
strongly correlated quantum fluids and could be used to 
investigate their properties. 



A quite different ID system where one may consider 
to experimentally realize a sonic black hole and to test 
Hawking radiation would be a system of charge parti- 
cles, like electrons in quantum wires or ions in elongated 
traps. The long-range part of the Coulomb interaction 
could substantially renormalize the Hawking tempera- 
ture in long wires in a way similar to how it does for 
the scattering from an impurity p^ . 

Conclusion — A Fermi-dcgcncratc flow of particles 
propagating in a ID channel in the presence of a smooth 
effective barrier provides an exact solvable quantum me- 
chanical model for the sonic analogue of Hawking radi- 
ation. The agreement with Hawking's calculation is ob- 
tained in the limit of a very smooth barrier, which is the 
equivalent in Hawking's derivation to a slowly evolving 
event horizon. The thermal and incoherent characteristic 
of Hawking radiation is discussed in terms of a measur- 
able quantity such as the dynamic structure factor. 
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